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A METHOD FOR DETERMINXNG FIELD RADT ATTQM L EVEI/S FOR A 

RADIATING DEVICq 
FIELD OF THE INVENTION 

The present invention relates to radiating 
devices such as antennas. 



BACKGROOm) OF THE INVENTION 



For health and safety reasons radiation levels 



from antennas must be accurately monitored and controlled. 
Regulatory bodies require radio frequency wireless 
10 operators to ensure radiation levels comply with the 
relevant health and safety standards in place. This 
requires the collection of technical data relating to 
radiation levels and exposure limits around a relevant 
antenna , 

15 Existing technig:ues for radiation estimation 

approximate the characteristics of the antenna to a point 
source and use far field theory to calculate radiation 
levels for an antenna. 

The present invention is aimed at providing an 

20 iii^roved method for determining radiation levels for 
radiating devices such as antennas. 
SUMMARY OF THE INVENTION 

According to one aspect of the present invention 
there is provided a method for determining field radiation 

25 levels for a radiating device comprising the steps of 
determining far field radiation characteristics of a 
radiating device, providing a model of the radiating 
device, which model approximates the determined far field 
radiation characteristics and determining a near field 

30 radiation characteristic from the model for at least one 
point in space. 

It is preferred that the model is chosen to 
approximate near field radiation characteristics. 

It is preferred that the method includes the step 

35 of determining a boundary between near field and far field - 
radiation. 

The method may include the step of determining 
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near field radiation density from the model. 

Preferably the method includes the step of 
determining near field radiation power density level over 
a plurality of positions in space. It is to be understood 
5 that space includes occupied (physical structure present) 
or unoccupied space and any particular area around the 

radiating device. 

It is preferred that the method includes the step 
of radiation pattern or gain characteristics of the 
10 radiating device from the two orthogonal far field 
radiation patterns. 

It is preferred that radiation pattern or gain 
characteristics are determined from documented data such 
as that available in handbooks. 

Preferably the method includes the step of 
determining 3dB beam width in the two orthogonal far field 

radiation patterns. 

It is preferred that the method includes 
determining physical characteristics of the device to 
20 determine the far field radiation characteristics. 

Preferably the method includes the step of 
providing a model including representing the device by a 
plurality of radiation sources. 

According to one embodiment of the invention the 
25 radiating device comprises a wire antenna. 

The method may include the step of providing a 
model including estimating the length and spacing of each 
wire element forming the wire antenna. 

Each radiation source may comprise one wire 

30 element. 

Preferably the method includes calculating mutual 
coupling between all wire elements of the radiating 
device . 

Preferably the method includes assembling an N by 
35 M impedance matrix for the radiation sources. 
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Tbe method preferably includes calculating 
voltage for each wire element. 

The method may also include determining current 
in each wire element • 
5 Preferably th^ method includes multiplying an 

inverse inipedance matrix by a column voltage vector (zero 
for parasitic elements and 1 volt for driven elements) . 

The method preferably includes assigning one or 
more Huygen's wavelet point sources to each wire element. 
20 The method may also include calculating magnitude 

and phase of each point source from ctirrent determined. 

The method preferably includes assxming a sin (9) 
dependence for each wire element ((9) being measured from 
the direction of the elements) and summing the 
15 contribution of each point source to each point in space 
within the near field radiation pattern* 

It is preferred that the method includes both 
near field and far field effects when calculating the 
contribution of each point source. 
20 The method preferably includes scaling field 

strengths determined at each point in space, by power 
supplied to the radiating device. 

The method preferably includes radiating devices 
which are Yagi-Uda^ log periodic, single or phased arrays 
25 of monopoles, dipoles, rhombic antennas and other regular 
or irregular wire antennas. 

It is preferred that the method includes 
providing a single point source for each wire element with 
a length less than half the wavelength of radiation 
30 emitted from the radiating device. 

The number of point sources for a radiating device will 
preferably be a multiple of the wavelength divided by two, 
i.e. twice the length of the radiating device divided by 
the wavelength. 

35 Preferably the method includes providing scaling 

factors for miscellaneous effects such as weather, 
obstacles, other antennas, metallic stanictures and 
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dielectric structxires or other factors which affect 
radiatioB characteristics. 

According to another embodiment of the present 
invention the radiating device is an aperture antenna. 
5 Preferably the method includes determining 

physical characteristics of the aperture antenna from the 
beam width characteristics. 

The physical characteristics include physical 
size, impedance, the size of the aperture and field 

10 distribution. 

The step of providing a model may include 
representing the aperture by at least one Huygen's wavelet 
soiirce . 

Preferably the method includes the step of 
15 summing the contribution from each wavelet source at each 

point in space. 

It is preferred that the method includes the step 
of summing the contribution from each wavelet source over 
a three dimensional coordinate system in space, e.g. 
20 rectangular, circular or polar coordinate system. 

Preferably the contribution from each wavelet 
source determined includes power, voltage and current. 

It is preferred that power density level at each 
point in space is determined using known power density 
25 formulas such as those described in more detail 

hereinafter . 

According to another aspect of the present 
invention there is provided a method of estimating 
radiation power density of electromagnetic radiation 

30 comprising the steps of identifying a radiating device, 

representing the radiation device as a plurality of point 
sources which radiate electromagnetic radiation, 
estimating power density level at a plurality of positions 
in space for each point source and determining the total 

35 power density level at each position, by summing the 
contribution of each point source to the respective 
positions in space. 
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It is preferred that the method includes 
displaying the power density level for a plurality of 
positions* 

Preferably the method includes displaying the 
5 power density levels in graphical, tabulated, 
diagrammatical, pictorial or other form. 

The method may include choosing positions which 
fit into a two-dimensional or three dimensional coordinate 
system. 

IQ Preferably point sources include any portion 

which requires at least two portions to represent the 

radiating device* 

Preferably the method includes summing the power 
density level determined at one position for all point 
15 sources representing the radiating device prior to 

estimating the power density level at another position. 

Alternatively the method includes estimating the 
power density level for the plurality of positions for one 
point source, storing the estimated power density levels 
20 then estimating and storing the power density levels for 
the plurality of positions for another point source and 
summing stored power density levels for each point source 
at each position to calculate a resultant power density 
level at each position. 
25 It is preferred that the method includes 

identifying a plurality of radiating devices. 

The method may include representing each 
radiating device as a plurality of point sources. 

Preferably the method includes calculating the 
30 total power density level for each radiating device. 

The method may include calculating the total 
power density level for each point source at each position 
and summing power density levels calculated at each 
position. 

35 Preferably the distance between each point source 

is determined by the distance between points in space. 
Thus in a rectangular coordinate system the spacing 



wo 03/10482S 



- 6 - 



PCT/AUOl/00700 



10 



15 



20 



between grid points will determine the spacing between 
point sources, with closer distances between adjacent 
points in space resulting in more point sources for the 

radiating device. 

It is preferred that the electromagnetic 
radiation measured is radio frequency electromagnetic 
radiation. However other radiation is included. 

Preferably the method includes calculating far 
field and near field tapering characteristics for the 
radiating device. This may be accomplished using 
predetermined formula or field measurements. 

Preferably the error measured is radio frequency 

Preferably the method includes calculating the 
far field and near field tapering characteristics for each 
position. 

Preferably far field Pd formula is 



GM.\5 

PoweratAntenna * 1 0 
47iDi^ 



Gd = Antenna gain with respect to dipole of 
analysis single 

PaA = Power sent to antenna after lossy items- 
Di = Distance from antenna. 
Preferably D far field distance = 

25 ~ if aperture > I ^ 

»tions - define break poxnt of 

.16 times ifmrture>n 

Pdp„aOT = 41.3*Pd (para - parabolic aperture) 
Breal^oint = -16 Dff - parabolic aperture 
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= .25 Dff - rectangular aperture 

According to smother embodiment of the present 
invention there is provided a method of determining field 
5 radiation levels for a radiating device comprising the 
steps of determining far field radiation characteristics 
of a radiating device, determining the boundary between 
near field and far field radiation, determining the 
displacement of a point in space relative to the closest 
10 point on the radiating device and calculating the power 
density level at the point in space. 

Preferably, the method of determining field 
radiation levels for a radiating device utilising the 
closest point method includes the step of determining 
15 characteristics of the radiating device in a similar 
maimer to that used in connection with any one of the 
previous aspects or embodiments of the invention. 

It is preferred that the power density level is 
determined by the aforementioned formula for Pd. 
20 It is preferred that the power density formulas 

are modified according to modification factors affecting 
the gain of the radiating device and the degradation of 
radiation as a function of the displacement of the point 
in space from the radiating device. 
25 The method may include modeling the radiating 

device as a plurality of point sources. 

Preferably displacement is determined by 
determining X, Y, Z vectors in space. 

Preferably the displacement is determined using 
30 azimuth and elevation angles of the point in space. 

It is preferred that the method includes 
determining the orientation of the radiating device in 
space to determine down tilt. 

The method may include the step of determining 
35 any offset in the displacement as a result of down tilt o 
the radiating device. 

The method preferably includes the step of 



wo 03/104825 



PCT/AU03y00700 



determining if the point in space is outside the width 
plane or length plane or height of the radiating device. 

Preferably the method includes the step of 
calculating the effective reduction of antenna aperture as 
5 a result of the displacement of the point in space from 
the radiating device. 
BRIEF DESCRIPTION OF THE DRAWINGS 

A preferred embodiment of the present invention 
will now be described by way of example only with 
10 reference to the accompanying drawings: 

Figure 1 shows a schematic diagram of a method of 
determining power density of a point in space determined 
using a closest point algorithm; 

Figure 2 shows a schematic diagram of a method of 
15 performing power density of a point in space determined 
using a Hyugen's wavelet method for an aperture antenna; 

Figure 3 shows espostire limit boundaries of a 
test antenna using traditional modeling techniques; 

Figure 4 shows exposure limit boundaries for the 
20 same test antenna using a closest point algorithm 
techni<2ue; 

Figure 5 shows a power density plot for the test 
antenna exposure limit botmdaries shown in Figure 3; and 

Figure 6 shows a power density plot for the test 
25 antenna modeled according to the closest point algorithm 
in accordance with the present invention. 

The best mode of performing the present invention 
will be described with reference to three different 
embodiments. 

30 Antennas typically fall into two main categories. 

One of these is wire antennas and the other is aperture 
antennas . 

A method for determining power density radiation 
levels for aperture antennas in accordance with the 
35 preferred embodiment of the^ invention may incorporate a 
closest point algorithm technique or a multiple point 
source technique. 
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POWER DENSITY CALCULATIONS 

AS described in the Microwave Engineers Handbook, 
Volume 2, Artech House 1971 to Sadd, Theodore [1], the 
simplest way of estimating the power density radiated from 
antennas is to apply the far field power density formula 
to a point source representation of an antenna. 

To achieve sufficient accuracy manufacturers far 
field gain patterns must be used. These exhibit the 
antenna far field gain characteristics for all directions 
(i.e. 0 to 360") in the horizontal and vertical planes. 
The far field power density formula is given by 



„ . PoweratAnteima* 10 '* 
Pd= 7-:^ 



Where Pd is the estimated power density, Gd is the 
15 antennas gain with respect to a dipole at the analysis 
angle, PoweratAntenna is the power sent to the antenna 
after lossy items such as a signal feeders, Di is the 
distance from the antenna. Units for the formula are 
Watts per centimeter s(iuared. 
20 The inventor has observed that this far field 

estimation grossly overestimates the power density in the 
near field area around the anteima. From observations and 
theory, according to one embodiment of the invention an 
algorithm has been developed for determining how far away 
25 the far field is from an antenna (and hence when the far 
field power density formula becomes accurate) . The 
algorithm is described in greater detail hereinafter. 

If the effective aperture of the antenna is less 
than 1 wavelength at the operating frequency then: 



wo 03/104825 PCT/AU03/0D7O0 

- 10 - 

wavelength 

If the effective aperture is greater than or 
equal to 1 wavelength at the operating frequency, then: 



2*Aea 
5 ^^.r-M' - wavelength 

wherein Aea is the emtenna's effective apertxxre. 

NEAR FIELD CALCUIATIONS 

To coinpensate for the overestimating effect of 
just using far field calculations it is necessary to 

10 deterittine the reduction in gain in the near field due to 
the finite size of the aperture. Prom power density 
observations for uniform line sources and tapered 
illumination aperture antennas a general break point 
distance must be determined where calculations change 

15 modes between near field and far field. An associated 

taper method can be used to calculate near field radiation 
levels. The break point for parabolic antennas is 
determined as 0.16 times the far field distance. The 
taper method for the parabolic antenna is defined as: 

20 W^,™,„,«.,=4).3*Prf 

where the Dl used in the Pd calculations is just Pfar- 
field, irrespective of how close a distance to the 
antenna. The break point for rectangular aperture 
antennas is defined as 0.25 times the far field distance. 

25 The taper method for the rectangular antenna is defined 
as; 

where the Di' used in Pd is replaced with Di times the 
breakpoint distance. This effectively makes the power 
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density decay with a rate of 1/Di in the near field rather 
them (which occurs in the far field) . 

These taper methods can be seen to still 
overestimate in the near field, but to a much lesser 
5 extent. This small overestimation builds some tolerance 
into the calculations to cater for anteimas that may have 
differing properties different to those exhibited, 
APERTURE ANTENNAS 

Before power density levels can be determined for 
10 an antenna it is necessary to identify key characteristics 
of the antenna. This enables the antenna to be accurately 
modeled* 

According to one embodiment the far field 
radiation pattern can be determined from the Fourier 

15 transform of the electric field across the aperture. When 
most of the radiation is concentrated in a highly 
directional beam, then close to the antenna there is no 
significant diffraction pattern when observed close to the 
main beam direction. That is radiation from any point 

20 source in the aperture travels almost the same distance as 
radiation from any other point source in the aperture when 
the observation point is confined to the non- shadow 
region- the aperture dimensions translated in the 
direction of the radiation, a sixiiple inverse square law 

25 can be applied and a good approximation to the radiation 
levels is obtained. In the shadow, i.e. outside the 
aperture translated in the direction of the main beam, 
diffraction effects occur and the field is attenuated. 

The optical equivalent is a parallel ray of light 

30 falling on an irregular hole in an opaque material. When 
observed very close to the opaque material, the shape of . 
the light beam is the same as the irregular hole (the 
Fresnel region) . When observed far from the irregular 
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hole, the image is that of a circle (Fraunhofer region). 

In order to calculate the near and far field 
radiation density given only two orthogonal radiation 
patterns and the power supplied to the aperture antenna, 
5 the mechanical structure must be determined from these 
radiation patterns. Once an approximate aperture 
distribution is known the antenna is modeled by point 
radiators/ (^^Huygen's wavelet sources) across the 
aperture • 

IQ According to the closest point algorithm method a 

full three-dimensional analysis of the antenna is achieved 
over the antennas length and width • The closest point 
algorithm is set out as follows: 

1. Find the displacement of the analysis point 
25 and the centre of the radiating element. 

2. Calculate the effective reduction of 
antenna aperture as a result of the 
displacement. 

3. Find the displacement vectors y and z 
20 dimension values between the centre of the 

antenna and the point of analysis. 

4. Rotate the displacement vector by the 
antenna's mechanical downtilt in three- 
dimensional space (to calculate the 

25 analysis with respect to the antenna face) . 

5. Find the Azimuth angle (the angle in the XY 
plane from the centre of the antenna's 
aperture face to the analysis point in the 
X-Y plane); 

30 6. Find the Elevation angle (the angle in the 

Z plane from the centre o£ the antenna's 
aperture face to the analysis point in the 
Z plane ) ; 
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analysis height going along the antenna 
height plane (obtain through the rotated 
displacement's z dimension); 

8. Prom the rotated displacement, obtain the 
5 analysis distance along the antenna width 

plane (obtained through the rotated 
displacement's xy plane vector's component 
along the antenna face 

9. If the offset (analysis displacement height 
20 + electrical downtilt compensation is 

inside +/-1/2 the effective antenna height, 
then the new source Z position is set to 
the antenna's centre height + the Z 
component of the downtilted offset 
15 Else 

The new Z position is set to the sign of 
the offset *l/2 the effective antenna 
height; 

10. If the antenna width plane offset is inside 
20 =/-l/2 the effective antenna width, then 

the new source X and Y positions are set to 
the antenna's centre X and Y positions 
the width offset along the suxtenna face + X 
and y components of the downtilted 
25 elevation offset. 

Else, 

•The new X and Y positions are set to the 
sign of the offset *l/2 the effective 
antenna width. 

30 The new antenna source positions are then used in 

the calculation at this analysis point. 

In the summary the closest point algorithm 
determines that for any point in space either within or 
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outside the near field/far field radiation boundary, the 
power density level at the point in space is measured 
relative to only one point on the aperture antenna. This 
point is the closest point according to the displacement 
5 data which is determined from measuring the X, Y, Z 

coordinates of the point in space as well as the elevation 
angle and azimuth angle. For each point in space the 
radiating point of the aperture is merely shifted to the 
most appropriate position to become the closest point on 

10 the antenna to the point in space. 

Referring to Figure 1 for points in space lying in 
the Par field 10 power density calculations can be made by 
singly using the previously highlight power density 
formula. However for points in the Near field 12 inside 

15 the boundary 11 the taper f\mction must be applied to 
accurately determine power density level. This taper 
function linearly reduces the effective aperture size in 
length and width as the distance f rcaii the antenna 
increases, until the Par field distance is reached, at 

20 ^ich point the closest point algorithm merges back into a 
point source analysis. The Near field calculations are 
based on the modifying factors incorporated into the power 
density Near field formulas hereinbefore described. 

The closest point algorithm also takes account of 

25 the orientation of the antenna including its tilt. 

For any point 13 in space which is outside the 
face of the antenna aperture, the closest point source 
would be on the perimeter of the aperture 14. 

AS an alternative to the closest point algorithm 

30 method power density levels can be determined using an 
accumulated point source method. 

Thus as shown in Figure 2 a point in space 15 has 
a power density level which is determined by the 
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accuimilated effect of point sources 16 on the antenna 
aperture 17 • 

As with the closest point algorithm method the 
aperture antenna must first be modeled from the two 
5 orthogonal far field radiation patterns. The size of the 
antenna aperture is then calculated in conjunction with 
data on the physical size of the structure. 

Once the size and shape of the aperture has been 
determined the aperture is then represented by a number of 
10 Huygen's wavelet sources. The number of wavelet sources 
may be determined by using a single point source for each 
element with a length less than half a wavelength. For 
elements longer than this the number of point sources 
increases proportionately with the number of half 
15 wavelengths in the element. 

As an exazople the Yagi-Uda style antenna is 
discussed. The 3dB beam width in the two orthogonal 
radiation patterns are used to gain an estimate of the 
number of elements in the Yagi-Uda antenna. A standard 
20 configuration is then used to represent the antenna. 

•Phe sum of the contributions from each point 
source to every point in space in the vicinity of the 
antenna is used to determined the Near and Far field 
intensities. Thus as shown in Figure 2 the contributions 
25 from each of the point sources 16 are combined to provide 
a power density level at point 15. 

AS with the closest point algorithm the field 
strengths are scaled in accordance with the power supplied 
to the antenna. 
30 EXPOSURE LIMITS COMPARISON 

The following diagrams display the e3cposure limit 
boundaries {using 2 and 10 W/m^ for the red and yellow 
zones respectively) of the existing and new calculation 
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techniques for plan (horizontal) and elevation (vertical) 
boresignts of the test antenna, the Argue CTA610D-R. Note 
that in the traditional full modeling solution the antenna 
is identified as a point source located at its phase 
5 centre. See Figure 3 

POWER DENSITY COMPARISON 

The following diagrams represent the power density 
graphical representations of the existing and new 
algorithms for the test antennas, of the type Decibel 
10 DB580y (omnidirectional) . Note that in the traditional 
full modeling solution, the antenna is identified as a 
point source located at its phase centre. See Figure 4 

WIRE ANTENNAS 

The Huygen's wavelet method can also be applied to 

15 wire antennas. Such antennas can be classified as Yagi- 
Dda, Log Periodic, Phased Arrays of Monopoles and Dipoles 
and Rhombic. The fields are calculated using a single 
point source for each element with length less than half a 
wavelength. For elements longer than this the number of 

20 point sources increases proportionately with the nuaaber of 
half wavelengths in the elements. 

As with the aperture antennas before the wire 
antennas can be properly modeled it is necessary to obtain 
the 3dB bandwidth in the two orthogonal radiation 

25 patterns. Once tbis has been done an estimate may be made 
of the number of elements N in the antenna. 

With the niunber of elements equaling: 

N = — where 1 is the length of the antenna and A, is 

A 

wavelength 

30 the length and space of each of the wire elements can be 
determined • 

The method then requires a calculation of the 



wo 03/104825 PCT/AU03/00700 

- 17 - 

mutual coupling between all elements in the Eirray and the 
asseznbly of the N x N impedance matrix. 

To obtain tiHa current I in each element the 
inverse impedance Z matrix is multiplied by the column 
5 voltage V vector (0 for parasitic elements, 1 volt for the 
driven element) . 
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10 

This gives an example for determining the current 
for an n element array. 

A Huygen's wavelet point source is then assigned 

15 to the current locations of each element in the antenna. 
"The magnitude and phase of the point source is directly 
proportional to the current calculated previously. 

Assuming a sin (8) dependence for each wire element 
(9 is measured from the direction of the elements), the 

20 contribution of each point source is added to every point 
in space in the vicinity of the antenna. This includes 
both near and far field affects. Thus the effect of each 
point source on a point in space is calculated in 
determining the overall power density level at the point 

25 in space. This is similar to the approach used in the 
Huygen's wavelet method used for the aperture antezma 
algorithm. 

Field strengths are also scaled dependent upon the 
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power supply to the antenna. 

The effect of metallic and dieletric support 
structures can be modeled by point sources using image 
theory. These images contribute to the array of point 
5 sources in the model and the field intensity calculations 
are made as before. 

In its preferred form the above described 
embodiments of the invention are implemented by computer 
software. The software is preferably configured to store 
10 data relating to different types of antenna. This data 
would include manufacturers antenna pattern files in a 
number of different standard formats. In this way it is 
possible to compare an antenna being modeled with existing 
data. 

15 If an antenna being modeled does not fit an 

existing manufacturers antenna pattern file the software 
is able to receive measured data relating to antenna 
patterns and create an appropriate file. This can be 
added to a main database for storing antenna pattern 

20 files. 

Typical antenna pattern properties which are 
stored include pattern type, frequency for that pattern, 
system loss, resolution, linear averaged, pattern cut, 
pattern type, electric tilt and effective gain. 

25 In addition to the above files storing basic 

antenna patterns may include data relating to the 
horizontal beam width, the vertical beam width and the 
front to back ratio. It is preferred that once the main 
characteristics of an antenna have been characterised the 

30 antenna can be modeled according to one of the above 

described methods. A power density graphical or visual 
representation then can be obtained and the software is 
able to highlight the near field and far field radiation 
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patterns using different colours. These differences are 
highlighted in black and white in Figures 3 to 6. 

It is preferred that the sof twcure is able to 
select a number of resolution options so that the 
5 resolution of the power density patterns can be varied. 

It is also preferred that the software is able t 
record safe power density levels and restrict the images 
produced for near field and far field radiation to the 
safety values. This diagrams can be produced which show 
10 only the hazardous levels of radiation. 



